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1.  INTRODUCTION 


1.1  Objective 

The  long-term  objective  of  this  research  is  to  device  the  technology  of  trans-rectal  fluorescence 
diffuse  optical  tomography  (FDOT)  to  utilize  zinc-specific  fluorescent  marker  for  detection  of 
prostate  cancer.  The  specific  aims  of  this  research  are  centered  on  developing  methods  to 
overcome  the  unique  challenge  of  the  proposed  technology. 

1.2  Background 

Zinc  is  a  well-established  metabolic  marker  of  prostate  cancer  [1],  Benign  prostate  tissues  secrete 
zinc  in  the  form  of  Zn3Citrate2,  and  zinc  has  a  concentration  of  approximately  lOmM  in  prostatic 
fluid.  In  cancerous  prostate  tissue  the  zinc  secretion  ceases,  resulting  in  five  to  ten  folds  of 
reduced  level  of  zinc  in  the  associated  prostatic  fluid.  If  a  zinc-tagging  fluorophore  can  be 
administered  to  the  prostate,  the  concentration  distribution  of  the  fluophore  will  be  proportional 
to  that  of  the  free  zinc,  and  the  uptake  of  the  zinc-specific  fluorophore  by  cancerous  prostatic 
tissue  will  be  ideally  none  or  practically  much  less  than  that  by  normal  or  benign  prostatic 
parenchyma,  which  we  call  it  the  “reverse-uptake'’'’  case.  The  detection  of  prostate  cancer  by  use 
of  zinc-specific  fluorophore  then  needs  to  recover  a  weak  target  (or  targets)  of  interest  within  a 
strong  background  fluorescence,  which  is  contrary  to  the  convention  of  reconstructing  a  strong 
target  (or  targets)  of  interest  within  a  medium  of  weak-fluorescence. 

1.3  Specific  aims 

This  research  has  the  following  aims  to  achieve:  (A)  to  study,  algorithms  that  render  FDOT 
reconstruction  of  “reverse-uptake"  target  fluorescence  in  trans-rectal  imaging  geometry;,  (B)  to 
develop  methods  of  extracting  the  spatial  information  from  imagery  of  sagittal  trans-rectal 
ultrasound  for  prior-guided  trans-rectal  FDOT  reconstruction;  (C)  to  configure  trans-rectal 
ultrasound-coupled  FDOT_  detection  system  for  timely  acquisition  of  fluorescence. 

More  specifically,  trans-rectal  FDOT  of  zinc-specific  fluorescence  for  prostate  cancer 
detection  is  subjected  to  the  following  technological  challenges:  (1)  to  reconstruct  a  fluorescence 
distribution  wherein  the  target  has  a  “ negative  contrast’’  with  respect  to  a  strong  background;  (2) 
to  recover  a  fluorescence  distribution  in  an  endo-rectal  imaging  geometry  that  poses  significant 
spatiallv-variant  ( depth-dependent)  sensitivity  in  the  measurement.  As  the  experimental  FDOT 
study  is  to  base  upon  a  continuous-wave  system  that  we  developed  previously  for  DOT  imaging, 
an  additional  challenge  to  FDOT  exists,  which  is  (3)  to  rapidly  differentiate  fluorescence  signals 
excited  by  multiple  steady-state  sources .  Besides  the  above  three  FDOT  challenges,  another 
demanding  issue  is  to  export  the  spatial  profile  of  prostate  from  2-dimensional  imagery  of 
sagittal  trans-rectal  ultrasound^  (TRUS)  (versus  the  conventional  imagery  of  transverse  TRUS) 
for  the  spatial  prior  to  trans-rectal  FDOT.  These  needs  in  trans-rectal  FDOT  of  zinc-specific 
fluorescence  provoke  non-conventional  approaches  for  system  instrumentation  and 
reconstruction  algorithm. 

1.4  Overview  of  the  key  progresses  in  year  1 

During  the  first  12-months  of  this  study,  we  developed  novel  instrumentation  and  algorithm 
methods  for  overcoming  the  challenges  (2)  and  (3)  described  above,  and  studied  the  challenge  (1) 
described  above.  The  novel  instrumentation  approach  is  to  utilize  a  wavelength-swept  light 
source  for  swept-spectral-encoded,  sequential  source  delivery  for  rapid  DOT  that  is  to  be 
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extended  to  rapid  FDOT  measurements.  We  have  also  implemented  a  new  algorithm  of 
“geometric-differential-sensitivity”  framework  that  is  different  from  the  methods  previously 
demonstrated  [2-6]  for  mitigating  the  spatially-variant  sensitivity.  This  algorithm  has  been  tested 
on  DOT  and  is  being  extended  to  FDOT.  As  TRUS  is  to  be  integrated  to  FDOT,  region-wise 
reconstruction  of  ‘  ‘negati ve-contrast  by  use  of  the  spatial  prior  assumed  from  TRUS  is  also 
investigated  by  simulation. 

During  the  first  12-months  of  this  study,  the  PI  has  also  been  trained  by  taking  courses, 
presenting  in  conferences,  and  publishing  in  journals,  for  his  future  academic  career. 


2.  BODY 

2.1  Swept-spectral-encoded  sequential  source  excitation  for  DOT  and  FDOT 

The  imaging  system  for  this  FDOT  research  has  to  capitalize  the  trans-rectal  DOT  system  we 
developed  that  is  configured  for  steady-state  (DC)  detection  of  multiple  source-detector  pairs  by 
a  spectrometer  and  a  CCD.  Since  our  system  was  in  DC,  there  has  been  two  ways  to 
differentiate  multiple  source-channels  for  DOT  detection  as  indicated  in  Fig.  1:  (1)  a  time¬ 
multiplexing  method  by  mechanical  switching,  and  (2)  a  spectral-encoding  technique  by  parallel 
use  of  discrete  laser  sources  or  by  using  a  narrow-band  light  source  after  splitting  the  light  with  a 
grating.  The  method  (1)  can  be  applied  to  FDOT  detection,  however,  the  mechanical  switching  is 
limited  by  the  date  acquisition  the  speed;  the  method  (2)  can  reach  rapid  data  acquisition, 
however,  it  cannot  be  applied  directly  to  FDOT  measurement  [8-10]  because  the  fluorescence 
signals  excited  by  simultaneously  applied  DC  source-channels  of  slightly  different  wavelengths 
cannot  be  differentiated. 
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Works  for  DOT 

1 

Does  not  work  for  FDOT 

Aj  A,  A,  A,  A.  A. 


Wavelength  encoding 
—by  a  narrow-band  source 


Potentially  fast 


Fig  1.  Source-endoding  technique.  (1)  Mechanical  swithing  works  for  both  DOT  and 
FDOT.  (2)  wavelength-encoding  works  for  DOT,  but  not  for  FDOT. 


There  are  the  need  of  utilizing  our  existing  DOT  system  for  FDOT,  and  the  necessity  of 
rapid  FDOT  data  acquisition  in  the  trans-rectal  imaging  mode.  In  this  research  we  have 
demonstrated  a  swept-spectral-encoding  source-encoding  technique  that  works  well  for  DOT  and 
has  the  potential  of  rapid  FDOT  data  acquisition  of  steady-state  fluorescence  (steady-state  in  the 
sense  of  neither  in  frequency-domain  nor  in  time-domain),  it  can  make  our  current  detection  unit 
adaptable  to  FDOT  measurement.  The  idea  is  to  time-multiplex  steady-state  sources  that  excite 
the  fluorescence,  by  sweeping  the  wavelengths  of  the  sources  in  a  band  of  excitation  spectrum — 
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a  “swept-spectral-encoding”  technique.  This  technique  is  an  evolution  of  the  spread-spectral- 
encoding  technique  we  demonstrated  for  rapid  and  endoscopic  DOT  by  using  a  broad-band  light 
source  [11-12],  By  sweeping  the  source  wavelength,  the  source  channels  will  be  excited 
sequentially,  as  shown  in  Fig.  2,  so  the  fluorescence  excited  by  different  “steady-state”  source 
(again,  neither  being  modulated  for  frequency-domain  detection  nor  short-pulsed  for  time- 
domain  detection)  can  be  differentiated.  This  method  has  the  potential  of  reaching  higher 
acquisition  speed  than  mechanically  switching  source-channels  for  steady-state  fluorescence 
detection  in  the  geometry  we  use  for  prostate  imaging. 


t  IB 


Time  encoding 
—by  a  swept-souroe 


Potentially  fast 


Fig  2.  Swept-spectral-endoding  of  the  source  allows  sequencing  the  source  excitation 
rapidlly  thereby  differenting  the  emission  channels  rapidly  in  FDOT  mode. 

A.  System  configuration 

The  swept-spectral-encoding  method  is  first  demonstrated  in  DOT.  As  is  shown  in  Fig.  3(a),  the 
light  from  a  swept-source  is  dispersed  by  a  spectrometer  (spectrometer  #1).  The  dispersed  light 
beam  thereby  linearly  scans  across  the  exit  window  of  the  spectrometer  #1  and  is  sequentially 
coupled  to  the  source  channels  linearly  arranged  at  the  output  plane  of  the  spectrometer  #1.  At 
the  detection  end,  the  detector  signals  from  each  source  channel  is  dispersed  by  another 
spectrometer  (#2)  and  focused  onto  the  CCD  sensor.  The  CCD  captures  one  image  for  each 
source  channel  sequentially  when  a  specific  channel  is  illuminated. 


(a)  System  schematics  (b)  Constructed  system 


Fig.  3  System  illustration  for  the  swept-spectral-encoding  based  DOT 
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Figure  3(b)  shows  the  actual  system  fabricated.  A  4mW  Superlum  BS-840-02  tunable 
semiconductor  light  source  is  controlled  by  Labview  to  sweep  in  the  range  of  838nm  to  858nm  at 
an  increment  of  0.05nm.  The  highest  sweeping  speed  is  lOOnm/ms.  The  source  light  is  lead  out 
by  a  single-mode  fiber  and  collimated  by  Thorlabs  C230TME-B  (focal  length  4.51mm), 
producing  a  0.08  mm-diameter  light  beam  cross  section  for  coupling  to  the  1mm  diameter 
source-channel  fibers  of  the  DOT  applicator.  The  detection  is  done  by  a  PIXIS  512  CCD  camera 
with  12. 3mmx  12.3mm  imaging  area  through  a  SpectroPro  2300i  spectrometer  with  300mm  focal 
length.  Source-detector  synchronization  is  achieved  by  National  Instruments’  DAQ  card  using 
LabView.  For  each  imaging  cycle,  the  PC  first  sends  command  to  tune  the  source  light  to  the 
predetermined  wavelengths  at  which  maximum  light  power  is  coupled.  Afterwards,  an  external 
trigger  signal  is  send  to  the  CCD  camera  to  capture  detection  signals.  Limited  by  the  source 
power,  a  minimum  exposure  time  of  170ms  and  an  extra  150ms  CCD  readout  time  is  required 
before  tuning  the  source  light  to  the  next  source  channel.  Therefore  for  each  channel,  omitting 
the  PC  commands  transferring  time,  approximately  320ms  is  required,  which  is  equivalent  to 
approximately  0.5  per  second  frame  rate.  However  such  acquisition  duration  can  be  reduced  if 
the  source  is  more  powerful  and  the  CCD  has  faster  data  transfer  rate. 

We  used  an  axial-imaging  DOT  applicator  of  20m  in  diameter  [7]  to  demonstrate  the 
principle.  As  is  shown  in  Fig.  4,  the  8  source  channels  and  8  detector  channels  are  fabricated 
evenly  interspersed  around  the  perimeter  of  the  20mm-diameter  probe  cross-section.  However, 
one  source  channel  (marked  in  Fig.  4(d))  is  discarded  because  of  its  significantly  low  coupling 
efficiency  due  to  fabrication  defect. 


(a)  Photography  (b)  Probe  inner  strcuture 


(c)  Source-Detector 
Geometry  -  Side  View 


(d)  Source-Detector 
Geometry  -  Axial  View 


Fig.4  A  circular  endoscopic  imaging  geometry 
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(a)  Experiment  setup  illustration  (b)  Solid  phantom  targets 

Fig.  5  Experiment  setup  and  materials 
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The  setup  for  experimental  validation  is  shown  in  Fig.  5(a).  The  probe  (applicator)  is 
submerged  in  a  10”x  10”x5”  tank  filled  with  intralipid.  The  inner  walls  of  the  tank  are  painted 
black  to  minimize  the  reflection.  The  solid  phantom  targets  to  be  imaged  are  fabricated  with 
black  plastic  with  assumed  infinite  absorption  contrast  to  the  background  intralipid,  and  material 
with  known  optical  properties  of  jxa=  0.0056  mm'1  and  /if =1.03  mm'1.  As  is  shown  in  Fig.  5(b), 
side  lengths  of  the  phantom  cube  ranges  from  5mm  to  15  mm.  The  centers  of  the  targets  are 
placed  at  the  imaging  plane  of  the  probe  and  the  horizontal  stages  are  used  to  control  the 
displacement  of  the  targets.  Further,  the  probe  is  rotated  with  respect  of  the  targets  to  examine 
the  azimuthal  sensitivity  of  the  system. 


B.  Imaging  single  target  of  different  sizes 

The  system  is  first  examined  by  targets  with  varied  sizes  as  shown  in  Fig.  6.  The  side  lengths  of 
the  target  are  5,  7.5,  10,  12.5,  and  15mm,  respectively.  Black  plastic  materials  are  reconstructed 
with  obviously  higher  optical  absorption  properties  as  shown  in  Fig.  6(b)  compared  to  the 
known-value  phantoms  as  shown  in  Fig.  6(c).  For  all  five  target  sizes,  the  targets  made  with 
black  plastic  are  recovered  at  the  same  azimuthal  location  although  the  recovered  volumes 
decrease  with  respect  to  those  of  the  actual  targets.  For  the  known-value  phantom  targets,  the 
absolute  absorption  values  are  close  to  the  actual  ones  when  the  cubic  volume  is  sufficiently 
large.  However,  targets  with  volumes  less  than  10x10x10  mm3  are  overwhelmed  by  background 
artifacts. 
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Fig.  6  Experiment  on  recovering  single  target  of  different  size,  (a)  Target  location  and 
size  illustration;  (b)  Reconstruction  results  on  black  plastic  targets;  (c)  Reconstruction 
results  on  known-value  phantom  targets;  (d)  Absorption  coefficients  plotted  along  the 
dash  lines  in  (b);  (e)  Absorption  coefficients  plotted  along  the  dash  lines  in  (c). 
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C.  Imaging  single  target  at  different  deaths 

The  system  is  also  examined  with  targets  at  varied  depth  as  shown  in  Fig.  7.  A  10x10x10  mm3 
cube  fabricated  from  black  plastic  is  imaged  at  side-to-probe  distances  from  0mm  to  15mm 
(equivalent  center  depths  of  5mm  to  20mm).  Recovered  absorption  distributions  are  shown  in 
Fig.  7(b)  and  the  absorption  values  are  again  plotted  along  the  probe  radial  direction.  It  is 
observed  that  1)  target  deeper  than  15mm  center-depth  is  difficult  to  recover;  2)  the  recovered 
volume  decrease  as  the  depth  increases;  and  3)  similar  to  experiment  shown  in  Fig.4,  all  of  the 
target  centers  are  recovered  closer  to  the  probe  at  the  depth  of  approximately  3mm  instead  of  the 
actual  varied  depth  up  to  15mm.  The  depth  recovery  inaccuracy  will  be  discussed  later  when 
introducing  a  new  reconstruction  algorithm  for  better  target  localization  in  depth. 


Fig.  7  Experiment  on  recovering  target  at  different  depths,  (a)  Target  location  and  size  (b) 
Reconstmction  results  on  10x1  Ox  10mm3  black  plastic  targets  (c)  Absorption  coefficients 
plot  along  the  dash  lines  in  (b) 


D.  Imaging  single  target  at  different  azimuthal  positions 

We  also  examined  the  azimuthal  uniformity  by  setting  a  target  3mm  away  from  the  optical  array 
at  varied  angular  position  as  shown  in  Fig.  8.  The  target  depths  are  desirably  recovered  in  all 
positions  studied  as  the  target  is  located  at  the  most  sensitive  region  of  the  imaging  geometry. 
Absorption  values  are  plotted  along  the  3mm  deep  circle  off  the  probe  (Fig.  8(c)),  as  is 
delineated  in  Fig.  8(b)  by  the  gold  dash  lines.  Although  the  values  fall  into  the  same  order,  a 
maximum  of  4  times  absolute  values  difference  is  observed.  The  azimuthal  non-uniformity  is 
mainly  related  to  imperfection  in  probe  fabrication. 


Fig.  8  Experiment  on  system  sensitivity  along  azimuthal  direction  (a)  Target  location  and 
size  illustration  (b)  Reconstmction  results  (c)  Absorption  coefficients  plot  along  the  dash 
lines  in  (b) 
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E.  Prospects  of  swept-spectral-encoded  sequential  source  excitation  for  FDOT 
We  have  successfully  demonstrated  the  swept-spectral-encoded  sequential  source  excitation  for 
DOT.  By  filtering  the  emission  light  from  the  excitation  band  at  detection,  the  method  is  readily 
applicable  to  FDOT,  which  is  to  be  demonstrated  in  the  future  stages  of  the  project. 

2.2  Geometric-differential-sensitivity  (GDS)  algorithm  to  improve  depth  localization  in 
endoscopic  imaging  geometry 

Due  to  the  unique  geometry,  endoscopic/trans-rectal  DOT  (and  hence  trans-rectal  FDOT)  is 
subjected  to  arguably  the  highest  level  of  sensitivity  variation  in  depth.  The  sensitivity  variance 
causes  miss-location  of  a  target  in  depth  towards  the  location  of  highest  sensitivity.  We  have 
demonstrated  a  new  algorithm  of  “geometric-differential-sensivity”  (GDS)  to  improve  the  depth- 
localization  in  endoscopic  DOT.  The  method  is  extendable  to  FDOT. 

A.  Principle  of  GDSJor  DOT  image  reconstruction 

The  principle  of  this  GDS  method  is  illustrated  conceptually  in  Fig.  9,  for  an  axial-translumenal 
geometry  with  8  sources  and  8  detectors  evenly  interspersed  along  the  array-medium  interface  as 
shown  in  (a).  In  this  geometry,  the  sensitivity  to  an  optical  heterogeneity  (e.g,  of  absorption) 
versus  the  depth  is  represented  by  the  dashed  curve  shown  in  (b),  with  a  near-boundary 
maximum  followed  by  nearly-exponential  decay  in  depth.  When  the  image  reconstruction  is 
based  upon  such  spatially-variant  sensitivity,  a  target  located  at  a  greater  depth  (hence  a  smaller 
sensitivity)  is  either  not  to  be  recovered  at  all  or  to  be  reconstructed  incorrectly  in  a  near¬ 
boundary  location.  Alternatively  for  an  optical  heterogeneity  at  a  fixed  depth,  the  shown 
sensitivity  profile  resembles  the  sensitivity  of  a  series  of  detectors  along  the  circumference,  with 
respect  to  the  same  source  that  is  azimuthally  closest  to  the  heterogeneity  of  interest.  A 
differentiation  between  the  sensitivities  of  two  neighboring  and  source-sharing  source-detector 
pairs  to  the  same  heterogeneity  carries  the  same  information  as  the  non-differentiated  sensitivity 
does,  however  with  the  variation  over  the  depth  significantly  reduced,  as  is  shown  by  the  solid 
curve  in  (b). 


-^-0%  Discarded 

,7  O  6  ’TT  _ 


Fig.  9  Conceptual  illustration  of  the  principle,  (a)  Imaging  geometry;  (b)  Magnified  area 
of  the  dotted  box  in  (a);  Depth-dependences  of  the  conventional;  (c)  and  geometric- 
differential;  (d)  sensitivity  profiles,  respectively;  (e)  Illustration  of  <i>:  '  in  equ.  X. 
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Analytically,  the  conventional  objective  function  to  be  minimized  during  DOT  image 
reconstruction  for  a  pair  of  source  i  and  detector  j,  denoted  by  {/,  /'},  can  be  expressed  as 


= 


0> 


hJi 


<D 


hj) 


J\x2 J.N  '  AVI 


(1) 


where  and  o;/'1  are  the  calculated  and  measured  photon  fluence  rates,  respectively,  SjU2, 


'Nxl 


is  the  updates  of  the  optical  properties  (assuming  absorption  coefficient  / ia  and  reduced  scattering 
coefficient  /A  only)  at  the  total  N  nodes,  and  J^n  is  the  conventional  sensitivity  matrix 
associated  with  source-detector  pair  {i,j}.  We  now  construct  a  surrogate  sensitivity  matrix, 


J, 


f i.i.t) 

U2N 


y  1x2  N 


- J , 


{i.J+k} 

U2N  j 


i=l:Nos,  k=\:Nod-\,  j=\:Nod-k 


(2) 


where  k  is  the  difference  between  the  serial  numbers  of  the  detectors,  “Aos”  is  the  number  of 
source  channels,  and  “NocT  is  the  number  of  detector  channels.  We  therefore  call  this  method  a 
geometric-differential  sensitivity  (GDS)  approach.  The  image  reconstruction  is  then  to  optimize 
a  surrogate  objective  function  [13],  as  expressed  by 


AT 


{'.Ml  _ 


{i.M)  ~ 

<D„  -<D„ 


A 


•  Aa 


(3) 


Note  that  in  equ  (2)  a  sign  change  of  the  resulted  term  j*'-  '•*<  ,  represented  by 
{'.  j,  k}  ^  -s  necessary  if  |k>Nod/2}  or  {k=  Nod/2  &  i<k  &  j<i}  or  {k=  Nod/2  &  i>k 
&  j+k>i}}.  The  profile  of  the  rearranged  surrogate  sensitivity  J,  after  column-wise  summation, 
is  like  the  one  shown  in  Fig.  9(d)  versus  the  profile  of  the  conventional  sensitivity  J  in  Fig.9(c). 

A  comparison  between  a  conventional  sensitivity  profile  and  the  GDS  profile  is  illustrated  in 
Fig.9(c)  and  (d),  calculated  using  NIRFAST  [14].  The  imaging  geometry  has  an  array  radius  of 
10mm  and  a  domain  radius  of  50mm.  The  finite  element  mesh  representing  the  imaging  domain 
includes  7708  evenly  distributed  nodes  and  15040  elements  with  the  Robin-type  boundary 
condition.  The  dynamic  range  of  the  introduced  GDS  profile  over  the  depth  is  approximately 
1/10  of  that  of  the  conventional  sensitivity  profile  over  the  same  depth.  The  nearly  depth- 
independence  of  the  surrogate  sensitivity  profile  is  expected  to  lead  to  much  improved  depth 
localization  when  compared  to  the  substantial  depth-degradation  of  the  conventional  sensitivity 
profile.  However,  due  to  the  increase  in  the  matrix  size  the  expected  improvement  in  depth- 
localization  by  GDS  approach  is  to  be  accompanied  by  an  increased  demand  in  computation  time. 


B.  Target  localization  by_  GPS-based  DOT  image  reconstruction — single  target  case 
GDS-based  endoscopic  DOT  image  reconstruction  is  examined  by  simulation  and  experiments 
for  two  sets  of  targets.  The  optical  parameters  employed  for  the  simulation  were  chosen  based  on 
experimental  studies.  Specifically,  the  background  had  jua= 0.0023  mm'1  and  //s-lmm'1 
corresponding  to  1%  bulk  intralipid  solution  at  840nm.  Optical  heterogeneity  was  introduced  by 
use  of  one  or  two  black-plastic  cubes,  whose  sharp  and  strong  contrast  over  the  background 
medium  was  necessary  to  evaluating  the  edge-spread-function  for  depth-localization.  The  cube- 
heterogeneity  was  treated  as  having  jua=  0.0115  mm'1,  or  five  times  of  absorption  contrast  over 
the  back-ground,  and  jus,=  1.1mm"1.  A  small  scattering  contrast  was  assigned  to  mitigate  the 
potential  cross-coupling  between  jua  and  jus\  A  1%  white  noise  was  added  to  all  forward  data. 

The  results  are  shown  in  Fig.  10.  The  array-facing-edge  of  the  cube  was  located  0,  5,  and 
10mm  away  from  the  optical  array,  or  equivalently  with  the  center  depths  of  6,  11,  and  16  mm, 
respectively.  The  results  of  the  conventional  and  GDS  methods  on  simulated  and  experimental 
data  are  displayed  column-wisely,  with  the  actual  location  of  the  cube  marked  with  dotted  square. 
In  addition,  the  radial  profiles  along  the  green  dotted  lines  crossing  the  cube  location  are 
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compared  for  each  row  in  logarithmic  scale,  as  is  shown  in  the  “Contour  Plots”.  It  is  observed  in 
simulation  that  both  conventional  and  GDS  methods  are  capable  of  recovering  targets  at  depths 
of  10mm.  However,  GDS  has  much  better  localization  of  the  depth  of  the  target.  In  the 
experimental  study,  the  recovering  of  the  target  at  edge-depths  of  0  and  5  mm  by  the  two 
methods  is  comparable  to  that  shown  by  simulations,  whereas  the  GDS  has  noticeably  better 
depth-localization  of  the  array- facing  edge  of  the  cube.  At  edge-depth  of  10mm,  neither  of  the 
methods  could  identify  the  target  at  the  correct  location,  however,  the  artifact  by  GDS  method  is 
negligible  comparing  to  that  by  the  conventional  method.  The  contour  plots  in  Fig.  10  imply  that 
in  both  simulation  and  experiment,  the  absorption  values  recovered  by  GDS  method  for  different 
target-depths  are  all  in  an  interval  between  e~6  to  e'5,  whereas  the  absorption  values  recovered  by 
the  conventional  method  vary  in  a  much  greater  range.  The  narrow  range  of  absorption  contrast 
demonstrates  an  additional  aspect  of  improvements  by  the  GDS  method. 


Fig.  10  Simulations  and  experiments  of  recoveimg  single  target  (Unit:  10’mm"1). 
Columns  (a)  and  (b):  simulation  results  by  conventional  and  GDS  methods, respectively. 
Columns  (c)  and  (d):  experiement  results  by  conventional  and  GDS  methods, 
respectively  .Blue  dotted  squares  indicate  target  positions.  Contour  plots  are  absorption 
profiles  along  the  green  dotted  lines  in  the  reconstruction  images.  Shaded  area  in  Contour 
Plots  indicates  the  target  locations.  GDS  improves  target  depth  localization  and  volume 
estimation. 


C.  Target  differentiation  by  GDS-based  DOT  image  reconstruction — dual-target  case 
Figure  1 1  demonstrates  simulation  and  experimental  studies  of  resolving  two  black-plastic  cubes 
of  10mm  edge-length.  The  array- facing  edges  of  the  cube  targets  are  placed  1mm  away  from  the 
probe  (corresponding  to  a  center  depth  of  6mm),  and  the  angular  separation  between  the  two 
targets  was  approximately  90  degrees.  The  medium  conditions  and  the  optical  properties 
assumed  for  the  two  cube-targets  are  the  same  as  those  in  Fig.  10.  The  actual  target  locations  are 
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marked  by  blue  dotted  line  in  Fig.  ll(a)-(d).  Along  the  azimuthal  cross  section  of  the  imaging 
geometry  as  is  marked  with  green  dotted  circles,  the  recovered  absorption  coefficient  profiles  are 
displayed  in  Fig.  1  l(a)-(d),  after  a  normalization  to  the  peak  of  the  extracted  profile  as  plotted  in 
Fig.  1 1(e).  Simulation  results  indicate  that  even  though  both  conventional  and  GDS  methods  are 
able  to  resolve  the  two  targets,  GDS  method  separates  the  targets  more  distinctly.  The 
conventional  method  based  on  the  experimental  data  failed  to  discern  the  targets,  yet  the  GDS 
method  resolves  the  two  targets  at  the  correct  locations  with  good  contrast  over  the  background. 


Fig.  1 1  Simluation  and  phantom  studies  on  resolving  two  targets.  (UnitlO'^mm'1).  (a)-(d) 
reconstruction  images;  Blue  dotted  square  indicate  the  target  positions;  (e)  Absorption 
contour  plotted  along  the  green  dotted  lines  in  (a)-(d),  shaded  area  indicating  the  targets. 

D.  Prospects  of  GDSzbased  FDOT  image  reconstruction 

We  have  successfully  demonstrated  a  method  of  mitigating  the  spatial-variance  of  sensitivity  in 
endoscopic  DOT  image  reconstruction  to  improve  target  localization.  We  are  conducting 
simulation  studies  to  compare  the  outcomes  of  this  new  algorithm  vis-a-vis  those  of  existing 
methods,  such  as  sensitivity-weighting  method.  Should  this  GDS  method  outperform  other 
existing  methods  in  target  localization  in  depth,  it  will  be  applied  to  endoscopic  FDOT,  as 
endoscopic  FDOT  is  subjected  to  the  same  level  of  spatial  sensitivity  variation  as  that  in 
endoscopic  DOT. 


2.3  Synthetic  study  of  “negative-contrast”  FDOT  with  spatial  prior 

Both  the  instrumentation  and  computation  modeling  works  in  section  1  will  be  extended  to 
fluorescence  imaging.  Currently,  we  have  investigated  in  simulation  region-wise  FDOT 
reconstruction  of  “negative-contrast”  targets  with  assumed  spatial  prior. 

A.  Region-wise  reconstruction  to  simultaneously  recover  emission-band  optical  properties  and 

the  fluorescence  yields  under  spatial  prior  for  negative  contrast  case 

Before  presenting  the  details  of  our  improvement  to  the  region-wise  reconstruction  algorithm,  we 
review  the  basis  of  FDOT  reconstruction.  Fluorescence  optical  tomography  usually  models  the 
propagation  of  excitation  and  emission  light  by 

V •  Dx (F)V4y  (F, ®)  +  [Max (F)  +  — ]4\ (F ,co)  =  (F, co)  (4) 

c 

V  •  Dm  (F)VO„  (F,  a>)  +  [juam  (F)  +  — ]0,„  (F,  co)  =  q0_m  (F,  co)  (5) 

c 

v  •  Dm  (F)VO»  m  (F,  co) + [pam  (F)  +  —  ]Ora  (F,  co)  =  0,„  (F,  co)qpaf  (F)  1  ^ 

c  l+[coT(r)Y 


10 


Where  subscript  “x”  and  “m”  represents  excitation  and  emission,  respectively.  ®  is  the  photon 
fluence  rate.  D  is  the  diffusion  coefficient  and  //a  is  the  absorption  coefficient.  /4/'  =£/rC//  is  the 
fluorophore  absorption  coefficient,  where  8/7  and  cji  are  the  extinction  molar  coefficient  and 
concentration  of  the  fluorophore,  respectively,  rj  and  r  are  the  fluorophore  quantum  efficiency 
and  life  time,  respectively,  qo  is  the  source  term.  F  and  co  are  the  pixel  position  vector  and 
modulation  frequency  of  the  source  light. 

Equations  (4)  and  (5)  are  usually  written  together  in  most  literatures,  as  the  differences 
only  existed  in  the  subscripts  of  “x”  or  “m”.  However,  more  accurate  FDOT  reconstruction 
should  include  the  following  procedures: 

•  To  implement  source  light  at  the  excitation  band  (qo.x)  and  collect  diffused  light 
signal  at  the  excitation  band  (®x  at  tissue  boundaries)  and  recover  absorption  and 
diffusion  coefficient  distribution  (juax  and  Dx)  at  the  excitation  band. 

•  To  implement  source  light  at  the  emission  band  (qo.m)  and  collect  diffused  light  signal 
at  the  emission  band  (®m  at  tissue  boundaries)  and  recover  absorption  and  diffusion 
coefficient  distribution  (juam  and  Dm)  at  the  excitation  band. 

•  To  implement  source  light  at  the  excitation  band  (c/o.x)  and  collect  fluorescence 
emission  signal  at  the  emission  band  (<Dm  at  tissue  boundaries)  and  recover 
fluorescence  yield  and  life  time  distribution  (rjjuaf  and  f)  at  the  excitation  band. 

Note  that  the  emission  band  signals  in  the  second  step  are  generated  by  external  source 
light  at  the  emission  band,  whereas  those  in  step  three  are  generated  by  fluorophore.  The  first 
two  steps  can  be  replace  by  multi-spectral  DOT  measurements  and  reconstruction  for  more 
reliable  determination  of  //ax.m  and  Dxm  values.  However,  such  arrangement  requires  external 
light  sources  at  more  wavelengths. 

As  region-wise  reconstruction  in  FDOT  is  over-determined  at  most  cases,  the  unknown 
parameters  in  step  two  (juam,  Dm)  can  be  solved  simultaneously  with  those  in  step  three  0///af  and 
r).  Such  arrangement  in  image  reconstruction  actually  reduces  the  complexity  of  FDOT  system, 
as  the  external  source  light  at  the  emission  band  is  no  longer  necessary.  The  reconstruction  is 
therefore  reduced  to  two  steps: 

•  Implement  source  light  at  the  excitation  band  (qo.x)  and  collect  diffused  light  signal  at 
the  excitation  band  (®x  at  tissue  boundaries)  and  recover  absorption  and  diffusion 
coefficient  distribution  (juax  and  Dx)  at  the  excitation  band. 

•  Implement  source  light  at  the  excitation  band  (qo,x)  and  collect  fluorescence  signal  at 
the  emission  band  (®m  at  tissue  boundaries)  and  recover  absorption/diffusion 
coefficient  distribution  (juam  and  Dm)  and  fluorescence  yield/life  time  distribution  (////af 
and  r)  at  the  excitation  band. 


We  only  utilize  DC  fluorescence  measurements  for  reconstruction,  co=0  and  t=0  in  equ. 
(4)-(6).The  Jacobian  structure  of  the  first  step  is  constructed  following  the  conventional  method 
and  we  structured  the  sensitivity  matrix  (Jacobian)  in  the  second  step  following  equ.  (4): 

<91n/u  01n/n  Sln/jj  <91n/u 

d/C,  dD„,  dWaf  dr 


5  In  Il2 

d/C, 


(7) 


gin  /77 

d/C, 


gln/?7 

dz 
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And  the  sensitivity  matrices  are  substituted  into  Levemberg-Marquardt  method  for 
iteratively  solving  the  unknown  values. 

At+i  =xk-  [JT  (xk)J(xk)  +  jLJY'jr  (xk  )v(xk  )  (8) 

The  region-based  reconstruction  is  simulated  in  Fig.  12  for  the  case  of  the  target  haing 
positive  optical  contrast  over  the  background  but  negative-contrast  of  fluorescence-yield,  and  in 
Fig.  13  for  the  case  of  the  target  having  negative  optical  contrast  over  the  background  and 
negative-contrast  of  fluorescence-yield.  The  imaging  domain  is  a  80mm x 80mm x 80mm  cube. 
Substructures  are  segmented  according  to  the  presumed  spatial  prior.  The  optical  array  includes 
7  sources  and  7  detectors,  arranged  in  parallel  rows  at  the  bottom  of  the  imaging  geometry.  The 
(b)-(f)  of  Fig.  12  and  13  shows  the  optical  properties  assignment  of  the  two  cases  of  medium- 
target.  We  assume  the  optical  properties  are  accurately  determined  in  step  one  and  only  the 
second  step  is  simulated.  /uax=0.006,  /r’sx=l.  The  fluorescence  quantum  efficiency  is  set  to  0.1  for 
all  regions.  The  reconstruction  initial  values  are  set  as:  pam  =0.01,  p’ sm  =1  and  ////af  =0.0001. 
Table.l  and  2  compares  the  assigned  and  the  reconstructed  values  of .  Note  that  the  scattering 
coefficients  ( p  ’sx,m  =1/(3-  Dx,m)-  /4x,m  )  are  shown  instead  of  showing  Dx.m  • 


Fig.  12  Simulation  geometry  and  procedures  with  region-wise  reconstmction  algorithm 
for  positive-optical-contrast  and  negative-fluorescence-yield,  (a)  Simulation  geometry; 
(b)-(d)  pax  and  p’sx  distribution,  reconstruction  parameters  in  the  step  1;  (d)-(f)  pam 
Ip  sm/^af  distribution,  reconstmction  parameters  in  the  step  2. 


Fig.  13  Simulation  geometry  and  procedures  with  region-wise  reconstmction  algorithm 
negative-optical-contrast  and  negative-fluorescence-yield.  (a)  Simulation  geometry;  (b)- 
(d)  pax  and  p  ’sx  distribution,  reconstmction  parameters  in  the  step  1 ;  (d)-(f)  pam  Ip  W7/4f 
distribution,  reconstmction  parameters  in  the  step  2. 
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In  Table  1  and  2,  although  the  reconstructed  values  do  not  agree  exactly  with  the  set 
values,  the  relative  positive  or  negative  contrasts  among  the  regions  are  accurately  identified. 


Table  1  Comparison  between  the  set  values  and  the  reconstruction  results  given  in  Fig.  12. 
Note:  The  optical  properties  in  excitation  band  is  assumed  as  known,  only  the  unknown 
_ parameters  in  the  second  step  (Fig.  10  (d)-(f))  are  recovered  here. _ 


regions 

/^am  / 

mm1 

ft  sm  / 

mm"1 

/^am  / 

mm"1 

ft  sm  / 

mm"1 

/V 

mm 

Set  values 

peripheral  tissue 

0.006 

0.8 

0.0072 

0.88 

0.0018 

prostate 

0.02 

1.25 

0.0240 

1.375 

0.0200 

rectum 

0.008 

1 

0.0096 

1.1 

0.0018 

tumor 

0.05 

2.0 

0.06 

2.2 

0.0018 

Reconstructed 

values 

peripheral  tissue 

- 

- 

0.0070 

0.88 

0.0018 

prostate 

- 

- 

0.0240 

1.444 

0.0194 

rectum 

- 

- 

0.0095 

1.104 

0.0018 

tumor 

- 

- 

0.04 

1.465 

0.0015 

Table  2  Comparison  between  the  set  values  and  the  reconstruction  results  given  in  Fig.  13. 
Note:  The  optical  properties  in  excitation  band  is  assumed  as  known,  only  the  unknown 
_ parameters  in  the  second  step  (Fig.  10  (d)-(f))  are  recovered  here.  _ 


regions 

/^am  / 

mm"1 

ft  sm  / 

mm"1 

/^am  / 

mm"1 

ft  sm  ! 

mm"1 

/V 

mm'1 

Set  values 

peripheral  tissue 

0.006 

0.8 

0.0072 

0.88 

0.001 

prostate 

0.02 

1.25 

0.0240 

1.375 

0.002 

rectum 

0.008 

1 

0.0096 

1.1 

0.001 

tumor 

0.004 

1.3 

0.0048 

1.43 

0.001 

Reconstructed 

values 

peripheral  tissue 

- 

- 

0.0072 

0.88 

0.001 

prostate 

- 

- 

0.0259 

1.23 

0.0024 

rectum 

- 

- 

0.0094 

1.11 

0.0010 

tumor 

- 

- 

0.0027 

1.22 

0.0020 

2.4  Future  works 

Future  works  include  the  following  aspects  of  the  study:  (1)  to  demonstrate  swept-spectral- 
encoding  sequential  source  excitation  for  FDOT  data  acquisition;  (2)  to  extend  GDS  method  to 
FDOT  reconstruction  of  both  “positive”  and  “negative-contrast”  cases;  (3)  to  configure  phantom 
and  experimental  setup  to  study  negative-contrast  fluorescence  detection;  (4)  to  develop  method 
of  extracting  spatial  profile  of  prostate  from  2-dimensional  imagery  of  quasi-sagittal  TRIJS. 


3.  KEY  RESEARCH  ACCOMPLISHMENTS 

•  Demonstrated  a  new  method  of  source-encoding  in  DOT  that  is  extendable  to  FDOT  for 
rapid  differentiation  of  fluorescence  excited  by  DC  sources. 

•  Developed  a  new  reconstruction  algorithm  in  DOT  that  is  extendable  to  FDOT  for 
mitigating  the  significantly  depth-dependent  sensitivity  in  trans-rectal  imaging  geometry. 

•  Implemented  a  two-step  FDOT  reconstruction  algorithm  to  the  synthetic  study  of  region- 
wise  reconstruction  of  “negative-contrast”  fluorescence  anomaly. 
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(2011). 

Conference  abstracts  and  presentations 

•  Xu  G,  Piao  D,  “Challenges  of  and  new  configurations  toward  fluorescence  diffuse  optical 
tomography  of  zinc-specific  biomarker  for  prostate  cancer  detection”,  poster  presentation 
in  2nd  Innovative  Minds  in  Prostate  Cancer  Today  (IMPaCT)  Conference  for  research 
funded  by  the  Department  of  Defense  (DOD)  Prostate  Cancer  Research  Program  (PCRP), 
Mar.  09-11,  2011,  Orlando,  FL.  Paper  PC094694-1920. 

•  Xu  G,  Piao  D,  “Near-infrared  diffuse  optical  tomography  based  on  a  wavelength-swept 
light  source”  Research  Symposium  2011,  Oklahoma  State  University.  The  presentation 
was  selected  to  the  First-place  of  oral  presentations  in  Biomedical  Science. 

Course  works  and  research  compliance  training 

•  The  PI  has  taken  courses  such  as  Optical  System  Design  and  Laser  Spectroscopy  offered 
in  Physics  department  for  additional  course-work  training.  He  earned  grades  of  A  from 
both  courses. 


5.  CONCLUSIONS 

In  conclusion,  progress  in  both  instrumentation  and  computational  aspects  pertinent  to 
trans-rectal  FDOT  of  negative  fluorphore  uptake  have  been  investigated  during  the  first  12- 
months  of  this  project.  These  progresses  are  related  to  the  tasks  specified  in  the  Statement  of 
Work  as  being  summarized  in  the  following: 

Task  1:  Synthetic  study  on  FDOT  of  fluorescence-strong  targets  within  a  weak  background  and 
FDOT  of  fluorescence-weak  targets  within  a  strong,  background  (partially  completed ) 

In  fulfillment  of  task  one,  a  region-wise  reconstruction  algorithm  simplifying  the  measurement 
procedures  in  conventional  FDOT  has  been  developed  and  validated  by  synthetic  studies.  A 
geometric-differential-sensitivity  reconstruction  algorithm  for  DOT  is  also  derived  for  the  depth- 
dependent  sensitivity  distribution  in  the  endoscopic  remission  imaging  geometry.  The 
effectiveness  of  the  algorithm  in  DOT  is  confirmed  by  both  computational  simulation  and 
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phantom  experiments.  The  effectiveness  of  the  algorithm  will  be  examined  in  FDOT  by 
simulation  and  experimental  studies. 

Task  2:  Optimize  image  processing  techniques  to  extract  spatial  prior  from  TRUS  for  3-D  FEM 
modej  generation  (postponed). 

As  the  instrumentation  works  in  first  year  focus  on  the  optical  part,  ultrasound  imaging  modality 
is  not  yet  integrated  into  the  system.  The  proposed  work  for  ultrasound  image  processing  is  to  be 
fulfilled  along  with  the  completion  of  the  task  3. 

Task  3:  Experiment  system  development  and  validation  (partially  completed, )_. 

A  system  based  on  swept-spectral-encoding  of  the  source  channels  is  demonstrated  for  DOT. 
This  novel  configuration,  with  the  optical  filtering  implemented,  will  allow  rapid  FDOT  of 
fluorescence  excited  by  multi-channel  DC  sources.  Further  studies  will  be  directed  toward 
experimentally  demonstrating  that  targets  of  negative  fluorescence  contrast  with  respect  to  the 
background  can  be  recovered,  with  the  investigated  instrumentation  and  algorithm  methods.. 
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